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PREFACE 


It has become apparent during the course of instruction 
at the Naval Postgraduate School that the phenomena of heat 
transfer is of ever increasing importance in many fields of 
engineering., The attainment of higher thermal efficiencies 
ultimately becomes a problem in heat trunsfer. The concevt 
of the heat transfer coefficient has become an accepted me- 
thod for predicting the relative ease with which heat flow 
፪ሽ оесиг. 

In conversations with Professor E. 5. Drucker, of the 
U. S. Naval Postgraduate School, the authors were made aware 
of the lack of adequate data on transfer coefficients as they 
relate to wet steam. Such lack of data was considered suf- 
ficiently important to be investigated. 

The experimental work of this thesis wés conducted at 
the U. S. Naval Postgraduate School from March to June of 1953. 

The authors desire to express their appreciation to Pro- 
fessor E. E. Drucker for his extensive assistance in the pre- 
paration of this thesis, to the Research Division of the Bab- 
cock and Wilcox Company for their invaluable aid in the in- 
strumentation of the test section, and also to Chief Jones, 
USN, and Chief Wallace, USN, for their aid in the construction 


of the experimental set-up. 
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ср Specific heat at constant pressure, Btu/(1b)(deg F). 

qd Diameter, feet. 

СО Nass velocity, 1b/(hr)(8a ft of” cross section). 

h Coefficient of heat transfer between fluid and surface, 
Btu/(hr)(so ft)(deg F). 

k Coefficient of thermal conductivity, Btu/(hr)(ft)(deg F). 

L Heated length of tube, feet. 

M Mass flow rate, lb/hr. 

ць КАЦІЇ уівсовіичу  16/(юг) (25). 

P Total pressure, 1b/sq in. 

тео .1116.... 

@ Fluid density, lb/cu ft. 

q Rate of heat transfer, Btu/hr. 

г Radius, feet. 

p Temperature, degrees Fahrenheit. 

v Fluid velocity, ft/hr 

x Moisture content of steam, lb water/lb mixture, percent. 

Nu i, Nusselt Group, dimensioniess. 

Br =>, Prandtl Group, dimensionless. 

Re че , Reynolds, Group, dimensionless. 

Subscripts: 

1 inlet. w water. t “tuber а depth 

2 outlet. S steam. о outer. i inner 





TABLE CF SYMBELS AND ABBREVLATIUNS 


Area or heat transfer surface, souare feet. 


. A! AJ መ ያ በ . 





. 


T ete aon ib om» 
= лы . 4 a с т, 2 APT lv ü 
— 


іже б бяюо ІРО (1 ЗИМ a 














om = jal ‹ሓ -“ 1 а 5 -. ከ) (የሸ መሠ 


лв ва оп ци ме 
. ۰ 0 "۱۰۱۲۷ ۶ Ты мн a 
жые гр S rr 
ЖҰ ¿dt? " 
amd ug лив „е 
she OA 

ሰ ኣኢኣጄቋ 


t x 


„ нс кай 
mq y то тето тё 

Mie am ДЈ] 

= ere Come) oh mort ere т 
DII IT ሙዛ 





















е 
i Me o u | n = 
ае ғ”. шә En 
(dm "9 تا‎ 
ilii Фо)» тв © ж ۳ 
ми Y MU мс % 





. 


SUMMARY 


The object of the development and preparation of this 
thesis was to determine the local or point coefficients of 
heat transfer to wet steam at moderate pressures, and to note 
their manner of variation with respect to moisture content 
of the steam. To accomplish this, steam at a known condition 
was passed through a vertical steel tube of 0.50 inches in- 
side diameter which was heated electrically. 

Coefficients were determined under the following con- 
ditions: moisture content in the range from O to Ó percent, 
absolute pressures up to 215 psia, and a variation in mass 
velocity from 73,000 to 172,000 pounds per hour per souare 
foot. 

The effects of moisture content were found to be signifi- 
cant, causing a rapid increase in the coefficient with in- 
creasing percent moisture. The local coefficients as they 
are herein presented vary "on то BTU's/hr-Ft*-°F for satur- 
ated steam to values approaching 6500 BTU! s/hr-Ft^-OF fior 


steam with a moisture content of about 6%. 
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INTRODUCTION 


Steam power has been long established as the primary 
means of ship propulsion in the U. 5. Navy. In the advent 
of atomic power for naval vessels, it is entirely probable 
that steam will continue to be the means by which this power 
is transmitted to the main engines. 

The Navy's interest in the improvement of the thermal 
efficiencies of boilers has led to the utilization of higher 
temperatures and pressures. It is felt that additional im- 
provement may be brought about by critical re-evaluation of 
the current concept of boiler design as it relates to heat 
transfer. Consequently, it is important to know as much as 
possible about the mechanism of heat flow to the working 
medium as it passes through successive stages from satur- 
ated water to saturated steam. One aspect of this problem 
may be solved by experimentally determining the coefficients 
of heat transfer to steam with varying degrees of moisture 
content. 

In this thesis work, experimentation was conducted in 
an effort to determine the coefficients for steam with a 
quality of 96 to 100 percent. Although this is a fairly nar- 
row range of moisture content, the results did tend to indicate 
the trend with which the transfer coefficients varied with 
quality. Time was a limiting factor in preventing further in- 


vestigation of steam with higher moisture content. 


226. 





፣ ም” 












«| 4 - -— 4 qued ars не тен 

— e жы. om m mama de 

ене ме (4 Ы. 3) b e Ram co fr 7 و‎ 

жн u en DLL и O “s 

жы. а пи 4 

н р 1. ре то s ° up. и 

" о шо сно Ян a де ፡፡ 55 

_— w m n m ‹6= 9 

| i wa. Rando cm UI mq 
| 








мов о سس‎ жі палата міне мо Lon nu 

| eee ሠ E dne bm 6d 
дамою Wo o lo oon cma cm: comes ሐ 8 
— pak: waman AAA MAIS то 
| ——— i በ. “<... a... v ዌ* የዘቹ 
Ve en ee ee 
—© | mao imt жә 


Ii um. tote ed -‏ © سسا ك 

|! чик cen n Ж СХС АБО. 

> т! == мес ее = Ww 
ЕХЕ 
-= kjee =p m ጠው cam mà 3 

~ MM І ТСЕ A IA 

. Cap m ሥመ መ o е —À ` APTA 




















NE 


Heat transfer coefficients of this kind, that is, coef- 
ficients relating to forced convection, are usually corre- 
lated with respect to the local Nusselt, Reynolds, and Prandtl 
numbers. Such correlation, however, is invariably made for 
fluids which are of a single phase only. Obviously this was 
not the case in the present investigation. In addition, a 
thorough survey of all existing literature failed to yield 
any data on values of density, viscosity, or thermal con- 
ductivity of wet steam, thus preventing the computation of 
the above dimensionless quantities. Further research to de- 
termine some of these parameters would be extremely useful. 

This work may be considered as an extension of the in- 
vestigation of heat transfer in two phase flow. It repre- 
sents some of the initial data available for steam at relatively 


low flow rates and low heat flux densities. 
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CHAPTER І 
REVIEW OF LITERATURE ON TWO PHASE FLOW 


A survey of the broad field of heat transfer indicates 
that a large quantity of data has been collected on the heat 
flow to fluids inside tubes. Unfortunately very little data 
has been published with respect to gas-liquid mixtures. 

The lack of information on temperature and velocity 
distributions has made it necessary to correlate experimental 
data on heat transfer coefficients by means of empirical 
equations in terms of certain dimensionless groups. Perhaps 
the most highly utilized equation is the Nusselt relation for 
single phase flow in forced convection through a tube (4):--- 

№=0.024 Re. pr0.4 (1) 

It has been found by McAdams and others (4), (5) Егош 
data for vaporization inside tubes that the local film co- 
efficients, in the preheat section where the fluid is ini- 
tially heated, are usually somewhat higher than would be pre- 
dicted from the Nusselt relation. This result is logical 
from the standpoint that some vapor bubbles may have formed 
in the preheating section, with the possibility that they 
disrupt the film locally as they are formed on the tube wall. 
As they move out into the fluid they would condense thus re- 
leasing their heat of condensation. The combined effect is 
an increase in heat transfer coefficient. In a comparable 


manner, it might be expected that as far as the influence on 
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the film coefficient is concerned, the coefficient will in- 
crease when gas passes through a system simultaneously with 
a liquid. 

In 1949 Bergelin(2) observed that if a gas and a liquid 
pass simultaneously in upward flow through a tube, as the 
gas-liguid ratio is increased, the following three flow con- 
ditions could be distinguished: 

(a) bubble flow, when gas bubbles pass individually 
through the tube; 

(b) slug flow, when alternatively slugs of gas and 
liquid pass through the tube, and | 

(с) annular flow, when the liauid flows in an annulus 
along the tube wall, while the gas passes at a much higher 
velocity through the center of the tube. Verschoor and 
Stemerding(8) in their investigation of the transfer coeffi- 


Cients to air and water mixtures verified the observations 


of Bergelin. They found that as their air to water mixture 
was increased, the heat transfer coefficient increased until 
a maximum was reached that coincided with the transition from 
slug to annular flow. Prior to reaching this maximum there 
was a change in slope in their curve which corresponded to 
the change from bubble flow to slug flow. 1% is interest- 
ing to note that their heat transfer coefficients increased 
with an increase in flow density, however this was of second- 
ary importance compared to the influence of varying the air- 


water mixture. 
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Yoder and Dodge@9) measured film coefficients of Freon-12 
boiling in a vertical tube. Their results showed that pres- 
sure and flow rate had very little effect on the coefficients, 
but as the relative amount of vapor increased from 40 to 100 
percent, there was a very significant decrease in the heat 


transfer coefficients. 





CHAPTER IT 
"^UIPMENT 


The equinment was a flow system in which steam was taken 
from a main line and passed through a throttle valve, a cooler, 
the test section, a condenser, and to a weigh tank. The es- 
sembly is shown schematically in Figure 1. The steam was sup- 
plied from a Babcock and Wilcox FM boiler at a pressure of 
200 psig. After passing through a throttle valve the steam 
was made to flow through a cooler which was intended for use 
as a method of varying the ouality. It was never subse- 
quently used since time limited the investigation to a per- 
cent moisture not greater than 6 percent. On either end of 
the test section an Ellison combination throttling and separat- 
ing calorimeter was installed for quality determination. Prior 
to the condenser an orifice of 0.166 inches diameter was in- 
stalled to maintain the pressure in the test section. A by- 
pass was provided around the orifice to permit adjustment of 
the flow rate independently of pressure. 

The test section as seen in Figures 2 and 4 consisted of 
a vertical steel tube having an inside diameter of 0.50 inches 
and an outside diameter of 1.25 inches. Its overall length 
was approximately 272 inches, but with an active length of 
2l inches for purposes of heat transfer.  Hecesses were cut 
in either end of the test section to minimize axial heat loss. 
The thermocouples in the test section were made of No. 28 


gauge chromel-alumel wire, and were installed by the Research 


-7- 





Division of the Babcock and Wilcox Company. Four pairs of 
thermocouples were installed along the tube, one of each pair 
at the outer surface of the tube, the other at a known depth 
below the surface thermocouples. The leads of the thermo- 
couples were threaded through a milled slot and out a i" 

tube welded to the test section. The original thought was 
that the thermocouples as located would permit determination 
of the temperature gradient across the wall of the tube, thus 
the inside surface temperature could be found by extrapola- 
tion. 

At either end of the test section, thermocouple probes 
were projected into the fluid stream for measuring the tempera- 
ture of the steam. These probes were obtained from Leeds and 
Northrop, and manufactured of iron-constantan wire. 

The heating element consisted of four independent heaters, 
each containing 29 8 feet of No. 18 Nichrome V wire. The 
wire was wound around the tube in a single layer. Each heater 
was approximately 6 inches in length. In order to prevent 
shorting of the heater coils, the wire was insulated with glass 
sleeving, and a thin layer of mica, 0.0015 inches thick, was 
wrapped around the tube before the heaters were wound. Power 
supply to the heaters was controlled by four 115v-l5amp Power- 
stats. See Figure 3. The total power to each heater was 
measured by portable wattmeters. 

The insulation of the test section was accomplished by 


winding the tube, outside the heaters, with strips of 1/8" 


обе. 
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thick asbestos cloth. This was followed by applying a length 
of standard 13" magnesia pipe insulation. The entire section 
was then covered with an aluminum sheet casing. Provisions 
were made so that air, heated to approximately the same tem- 
perature as the outer surface of the insulation, could be 
supplied to the space between the insulation and casing, thus 
minimizing any heat loss due to convection, and reducing the 


conducted heat loss through the insulation. 
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CHAPTER 111 
OPERATING PROCEDURE 


There were four primary variables associated with the 
apparatus: heat flux, flow rate, quality, and pressure. It 
was found, however, during preliminary experimentation, that 
the minimum moisture content of the steam available for use 
was approximately 6%. Interest at this time was primarily 
centered on the region between O and 10 percent moisture, 
thus necessitating throttling the steam to lower pressures 
in order to obtain a given quality at the entrance of the test 
section. Such an arrangement meant that only the heat flux, 
pressure, and flow rate could be varied independently of one 
another. The only other limiting factor associated with the 
equipment was, that due to the design of the heaters and the 
power supply used, the maximum heat input was restricted to 
1 Ки рег heater, or a total ot 4 kw. 

The procedure was to select an approximate quality at 
the entrance of the test section by setting the throttle valve 
accordingly. The flow rate was either taken as that provided 
by the pressure across the orifice in the line to the con- 
denser or adjusted by means of the by pass around the orifice. 
Power was then supplied to the heaters and the system was al- 
lowed to come to a steady-state condition. 

The approach to steady-state condition was determined 
by noting the change in the millivolt reading of the thermo- 


couples contained in the test section. When they became 
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approximately constant, steady-state was considered to exist. 
It normally took from 1/2 to 3/4 hours before the system 
reached equilibrium. It wes never possible to keep tempera- 
tures constant over a period longer than a few minutes, even 
though heat input and flow rate appeared unchanged. Each 

run was made at constant flow rate, heat flux, pressure and 
quality of the inlet steam. 

For each run, temperatures, pressures, flow rate, and 
power input were recorded. Data points were taken at various 
settings of power inputs, flow rates, and pressures. Most 
of the runs were concluded at a power input such that the 
minimum change in ouality which occurred through the test 


section was 3%. 
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CHAPTER IV 
METHOD OF CALCULATION 


The fundamental law of heat transfer by convection may 
be stated by the mathematical expression: (4) 

q +A IT (2) 
The factor of proportionality h, defined in the equation, is 
called "the film coefficient of heat transfér." It is this 
coefficient which was computed. 

For each set of conditions imposed on the test section, 
the density of heat flux, the inside surface temperature of 
the tube, the JT, the quality of the steam, and the local 
coefficient of heat transfer were calculated. Each of these 
factors was evaluated at four positions corresponding to the 
locations of the thermocouples inside the tube wall. 

The power dissipated in the heating elements was ob- 
tained from wattmeter measurements. The heat flux density 
was then calculated from the power dissipated and the dimen- 
sions of the heat transfer surface. 

Early in the process of collecting data for the investi- 
gation, it was discovered that the thermocouple pairs in the 
test section were indicating unreasonable temperature gradi- 
ents througn the tube wall. It became necessary to calcu- 
late the inner surface temperature by applying to Бим дербі 
thermocouple a computed temperature drop corresponding to 


the heat input. The temperature drop was calculated using 


210 





the Fourier conduction equation for steady-state heat trans- 
mission: (4) 
9 2 - ۸ ۸ له‎ (3) 
ФА 
For a circular tube, this equation may be integrated to give: 
Я = 2AM KE LAZ 
117 “Жж 
Тре depth thermocouples were shielded from any direct influ- 


ence of the heaters, and thus were considered to be more ac- 
curate than those on the surface of the tube. 


The temperature of the steam was measured at both the 
entrance and exit of the test section, however, these measure- 
ments were consistently erratic and failed to give desired 
results. It was felt that greater reliability would result 
from assuming the saturation temperature corresponding to the 
absolute pressure in the test section for computing the AT. 

In the preliminary phase of the investigation it was 
discovered that the reproducibility of the calorimeters, for 
determination of quality below the level which could be found 
by the throttling principle alone, was rather poor. With re- 
sults taken over a long period of time it was found that the 
average quality of the steam at the test section, at a pres- 
sure of 200 psig, was 94%. Further use of the calorimeters 
was discontinued, and this value of 94% for quality was used 
as a basis for all further computations. The quality of the 
steam, at each point in the test section corresponding to the 


location of the thermocouples, was computed by adjusting the 
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enthalpy of the incoming steam in an amount equal to the 
heat input to the test section. The quality at pressures 
lower than 200 psig was determined by the use of the Mollier 


Diagram. 
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CHAPTER V 
CONCLUSIONS AND KSSULTS 


One hundred runs were made of which forty-six were suc- 
cessful. The remainder of the runs were eliminated because 
the tube wall temperatures obtained indicated that steady- 
state conditions had not yet been reached. The ranges of 
pressure, temperature, quality, and mass rate of flow are 
given in Table I of Appendix І. 

In the process of collecting data a number of experi- 
mental errors became apparent. These errors, however, did 
not exceed the range of expectancy ordinarily associated with 
this type of investigation. The temperature determinations 
were expected to be in error at least 2 to 3%. This was not 
serious as long as the temperature difference between the 
fluid and the wall surface was large. At small temperature 
differences care had to be taken to insure that unreasonable 
results were not obtained. Another error is introduced in 
the calculation of the transfer coefficient by neglecting 
axial heat loss from the test section. In estimating this 
loss it was found to be small, in the range of 10 to 40 BTU/ 
hr. The error due to the radial heat loss was approximately 
2% of the input power. 

Curves for the heat transfer coefficient versus the per- 
cent moisture were made for each mass velocity. In addition 
a plot was made of the transfer coefficient versus the tem- 


perature difference between the wall and the fluid, and one 
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of the transfer coefficient versus the mass rate of flow. 
These curves are seen in Figures 6 through 12. The physical 
properties of steam-water mixtures were not available, thus 
preventing any correlation of the transfer coefficient with 
any existing data. 

The curves depicting the heat transfer coefficient 
versus the percent moisture show a very gradual rise in the 
coefficient up to moisture content between one and two per- 
cent. Above this value the coefficient increased rapidly as 
the moisture increased from two to five percent. A possible 
explanation of this phenomenon is that at high qualities 
there are only small droplets of water in the fluid, but not 
enough to materially disrupt the boundary layer between the 
steam and the wall of the test section. As the amount of 
moisture in the steam increases, the water droplets become 
larger and begin to rupture the boundary layer and collect 
on the wall of the test section, with a consequent increase 
in the heat transfer coefficient. The sharp change in slope 
of the curve is the point at which the water begins to col- 
lect in droplets on the tube wall. With a further increase 
in moisture, there is an approach to annular flow where the 
boundary layer becomes completely liquid. At this point 
the transfer coefficient should reach a maximum and then 
start to decrease. Fully developed annular flow was never 


reached in this experiment. 
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Figure 11 shows the influence of mass flow rate on the 
heat transfer coefficient. It is surprising to note that the 
coefficient decreased as the flow rate increased to a value 
of about 200 lbs/hr. With flow rates above this value, the 
coefficient increased. Only one point on the curve was es- 
tablished beyond this transition region. However, the pre- 
dicted curve is extrapolated (dotted line) on the plot. The 
initial decrease in the coefficient was possibly due to the 
water droplets on the wall being flattened out as the flow 
rate increased and thereby forming a liquid boundary layer. 
As the mass flow rate increased beyond 200 lbs/hr the veloc- 
ity of the gas was apparently sufficient to start wiping the 
water droplets off the wall instead of just flattening them 
out. This is how the continued upward trend in the curve 
was predicted. 

Figure 12 indicates curves of the transfer coefficient 
versus the difference in temperature between the tube wall 
and the fluid. This shows the correlation between the heat 
transfer coefficient and temperature difference, for a given 
range in quality of about five percent. These curves appear 
to be in the region of unstable film boiling when compared to 
the boiling curve for free convecticn(), but the coefficients 
determined in the investigation are necessarily higher than 
those of the boiling curve due to the influence of forced 


convection. 
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The available literature and results of this investiga- 
tion permits the conclusion that a maximum heat transfer co- 
efficient exists with respect to the moisture content of the 
mixture. This point has not been within the range of the 
subject investigation. The location of this maximum repre- 
sents an optimum for future design of certain types of heat 
transfer equipment. | 

Additional research would be amply justified in order 


that the peak in this curve be definitely established. 
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TABLE III 
CALCULATED TEMPERATURE DROP FOR VaRIUUS FOWER INPUTS 


Power Input per Temperature Difference Temperature Dif- 
Coil in Kilowatts Between Surface Thermo- ference Between 
couple and Inside Surface Depth Couple and 
Inside Surface 


300 11.94 2.96 
350 13692 3 226 
lOO 15.91 2772 
450 17.90 1,3280 
500 19.92 4 „66 
2210) e LA 5 
600 23.9 5200 
650 25.9 6.06 
700 27.9 6.58 
750 29.8 u.00 
800 31.8 1.46 
850 33.8 7.92 
900 35.8 8.40 
950 QNS 8.86 
1000 39.8 9.32 
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